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Abstract 

| We present a supersymmetric see-saw S4 model giving rise to the most general 

^ neutrino mass matrix compatible with Tri-Bimaximal mixing. We adopt the S4 x Z5 

Q_i, flavour symmetry, broken by suitable vacuum expectation values of a small number 

I of flavon fields. We show that the vacuum alignment is a natural solution of the most 

£T) ' general superpotential allowed by the flavour symmetry, without introducing any 

soft breaking terms. In the charged lepton sector, mass hierarchies are controlled 
by the spontaneous breaking of the flavour symmetry caused by the vevs of one 
i-C 1 doublet and one triplet flavon fields instead of using the Proggatt-Nielsen U(l) 

^p" 1 . mechanism. The next to leading order corrections to both charged lepton mass 

Oh| matrix and flavon vevs generate corrections to the mixing angles as large as 0(X C ). 

Applied to the quark sector, the symmetry group S4 x Z5 can give a leading order 
Vckm proportional to the identity as well as a matrix with 0(1) coefficients in the 
, Cabibbo 2x2 submatrix. Higher order corrections produce non vanishing entries 

£>• \ in the other Vckm entries which are generically of 0(X C ). 
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1 Introduction 



The Tri-Bimaximal structure (TBM) [T] of the neutrino mixing matrix is remarkably in 
agreement with the experimental results in the neutrino sector [2J. Within l-a error, the 
values of the mixing angles can be approximated by their TBM values [3]: 

tan 2 # 23 = 1 tan 2 12 = - sin 13 = . (1) 

The simplified structure of the mixing matrix suggested the possibility to be explained 
using some discrete non-abelian groups, added to the Standard Model, and containing 
a triplet representations to fit the number of lepton families observed in Nature. The 
symmetry A 4 [H [Sj (also in the context of Grand Unified theories [5]) emerged as a 
natural candidate because it is the smallest discrete group with triplet representation and 
it is sufficiently manageable to be broken differently in the charged and neutral lepton 
sectors, a necessary condition if we want to get a mixing matrix different from the identity. 
In the context of see-saw models, it is interesting to observe that the TBM structure in 
Ai is generally associated with a well defined relations among the complex eigenvalues of 
the light neutrino mass matrix: 

-!- = -!—?-, (2 , 

777,3 777i 7772 

implying that most A A models are quite predictive because of the reduced number of 
independent parameters. It should be stressed, however, that the realisation of the TBM 
strongly relies on the choice of symmetry breaking pattern. In fact, in the neutrino sector 
the group v4 4 is usually broken into a subgroup generated by the matrices 




and 

G 



which, in the basis where charged leptons are diagonal, leave invariant the most general 
neutrino mass matrix diagonalized by TBM: 





miight = I y x + v y-v \ . (3) 



It turns out that the representations of A^ only contain G and the invariance under £/ M _ r 
arises accidentally, as a consequence of the specific field content of the model. From this 
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point of view, the group S4 [T|-[S] arises as a natural candidate as a flavour group for neu- 
trino mixing because one can find a suitable representation of it containing simultaneously 
the previous elements. It should also be noted that the extension of the A 4 symmetry 
from lepton to quarks seems to be complicated by the absence of doublet representations, 
whose use is suggested by the heaviness of the top quark, which S4 possesses instead. 
However, the problem of reproducing small mixing angles and strong mass hierarchy in 
the up-sector at the same time can only be partially alleviated by the bidimensional rep- 
resentations of S4 because one or more fine tunings between the relevant Yukawas are 
invoked to correctly reproduce some of the m up / rridown quark mass ratios. In this paper, 
we build a constrained see-saw 54-based model for fermion masses and mixing which, 
compared to models already existent in the literature, realises the most general neutrino 
mass matrix diagonalized by TBM at leading order (LO). This is obtained allowing the 
right-handed neutrinos to couple to singlet, doublet and triplet flavon fields. The light 
neutrino masses depend on six complex Yukawa parameters and the typical A4 sum rule 
of eq.([2]) does not hold, leaving the model less predictive but more manageable. The 
mass hierarchy among charged leptons is obtained breaking the 5*4 symmetry by the vevs 
of a doublet and triplet flavon fields, without invoking any Froggatt-Nielsen £7(1) sym- 
metry. The unwanted couplings are forbidden imposing an additional Z 5 symmetry to 
the model. The resulting S4 x Z 5 symmetry is minimal from the point of view of the 
flavour symmetry and field content. We extend the S4 x Z5 symmetry to the quark sector 
including the left-handed components into triplet representations (and not into doublets, 
as usually done); we show that, even using a rigid structure like that proposed in this 
paper, some of the relevant features of the quark sector, like a good leading order Vqkm 
and quark mass ratios, can still be accounted for. The paper is organized as follows: in 
SectJ21 we discuss the relevant feature of the S4 symmetry and the structure of the model, 
presenting leading order results on neutrino as well as charged lepton mass matrices; in 
SectJ3]we compute the next to leading order corrections (NLO) to the vacuum alignment 
and the relevant higher order operators, both responsible for deviations from TBM mix- 
ing; in Sect JH we discuss some phenomeno logical results obtained from our model with all 
Yakawas constrained to be 0(1) and we also show that the whole model allows for ac- 
ceptable leptogenesis parameters. Sectj5]is devoted to the quark sector whereas in Sectj6] 
we draw our conclusions. 

2 The structure of the model 

We introduce here the structure of the model which leads to TBM in first approximation. 
We recall that S4, the permutation group of 4 objects, can be generated by the two 
elements 5* and T obeying the relations (a "presentation" of the group): 

S 4 = T 3 = 1, ST 2 S = T. (4) 

The action of the generators S and T can be assigned as follows: 

(1234) ^ s (2341) 
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(1234) — (2314) 



and the 24 elements of the group, belonging to 5 conjugate classes, are: 
1 



c, 
c 5 



S 2 = (3412), TS 2 T 2 = (4321), S 2 TS 2 T 2 = (2143) 

T, T 2 = (3124), S 2 T = (1423), S 2 T 2 = (2431), STST 2 = (4132) 

STS = (4213), TS 2 = (4132), T 2 S 2 = (1342) 

ST 2 = (1243), T 2 S = (4231), TST = (1432) 

TSTS 2 = (3214), STS 2 = (1324), S 2 TS = (2134) 

S, TST 2 = (2413), ST = (3142), = (3421), S 3 = (4123), S 3 T 2 = (4312) 



The inequivalent irreducible representations of S4 are 1\, I2, 2 and 3. It is immediate to 
see that one- dimensional unitary representations are given by: 



h 
1 2 



S 
S 



1 T 
-1 T 



(5) 



while the two-dimensional unitary representation, in a basis where the element T is diag- 
onal, is given by: 



T 



uj 
u 2 



S 



1 

1 



(6) 



Finally, the three-dimensional unitary representation is as follows: 
3 X : T = I I , S = i 




-1 


2a; 


2a; 2 


2w 


2a; 2 


-1 


2a; 2 


-1 


2a; 



(7) 



where a; = e 271 "*/ 3 = (— 1 + v / 3)/2, whereas in the 32 representation the generator T is the 
same but S is the opposite. It is useful to remind the product rules between the group 
representations: 



li®£ = £ 

l a <g>l 2 = li, 1 2 ®2 = 2, 
2® 2 = li©l 2 ©2, 2®3i 
3i®3 2 = 1 2 ©2©3 2 ©3 2 , 
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2 09 «3j = 

3i©3 2 , 3, 



3i 



where the indeces i,j = 1,2, with i 7^ j and £ indicates any other representation. The 
Clebsch-Gordan coefficients in the basis presented above are reported in Appendix A. 

The general neutrino mass matrix of eq.(j3]) can be obtained in the framework of the 
see-saw mechanism, 



nought 



-m T D m D 



(9) 



3 



where both Majorana (tum) and Dirac (mp) mass matrices are needed. They are derived 
from the most general lagrangian invariant under S4 x Z 5 and containing fields in any of 
the S4 representations, which are singlets, doublets and triplets. The group S4 is broken 
by means of suitable vev's of Standard Model singlet fields (flavons), whose alignments 
have to guarantee the correct entries of the Majorana and Dirac mass matrices. Group 
theoretical considerations help in understanding the pattern of symmetry breaking needed 
to generate the wanted matrix; in the representation of App.A, the elements S 2 and 
(TSTS 2 ) leave invariant the mu ght of eq.(J3]), that is 

S 2 rnught S 2 = mught 
(TSTS 2 ) mught (TSTS 2 ) = m Ught . 

This means that the vevs of the flavon fields in the neutrino sector have to be invariant 
under the subgroup generated by them. This strong condition for the flavon alignment is 
realized, for the triplet representation, by the field configuration 

(<ps) = v s (1,1,1). (10) 

For the bidimensional representation, the same matrices should leave the vev of a doublet 
field invariant; if we choose 

<A)= VA (1,1) (11) 

these correspond to the first matrix of the groups in the doublet representation. 
Therefore, we have found a scalar field configuration which remains invariant under the 
action of the same matrices that leave mu g ht invariant. In addition to the previous fields, 
we also include in the model a singlet flavon field £, with a non-vanishing vev. To realize 
the classical see-saw mechanism, we need to introduce right-handed neutrinos, which we 
assume to transform as a triplet representation of S4. To avoid large fine-tuning from 
terms of the form iWV c z/ c , like those discussed in [5], we properly tune the v c charge 
under z£]; whatever this Z5 charge is, if we want the singlet, doublet and triplet flavons 
to contribute to the neutrino mass matrix their charges have to be the same. We also 
attribute a Z 5 charge to the h u higgs boson in order to avoid the Weinberg operator 
5 = £h u £h u at the leading order [IT] , where £ is a triplet field of the SU(2) Standard 
Model lepton doublets. This also forbids a leading order Dirac mass term of the form 
(y c £) h u , which is then generated at the 0(1/ A) through couplings with the same flavon 
fields £, A and ips- The Z5 charges for the matter fields as well as for flavons and driving 
fields (discussed later) are reported in TabJTJ where we used the symbol u = e 2m ^ to 
indicate the Z 5 unit-charge. 

The lagrangian in the neutrino sector is then as follows: 

W v = jv c £ h u (y Ul ip s + y U2 A + y V3 £) + v c v c (a £ + b <p s + c A) , (12) 
1 Similar considerations, but using two discrete groups, have been discussed in |10j . 
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Table 1: Transformation properties of leptons, electroweak Higgs doublets and flavons 
under S4 x Z 5 and U(1)r . 



where a, 6, c and y Vi are complex Yukawa couplings. The Dirac mass matrix is obtained 
from the first term in eq. (TT2T) and it is given by: 

2 y Vl v s + Vu 3 u y U2 v A - y Ul v s y U2 va ~ Dv x v s \ 
Vu 2 v a ~ Vvi v s 2 y Ul v s + y U2 v a Vu 3 u - y Vl v s (13) 
Vv 2 va - y Vl v s y V3 u - y Vl v s 2 y Vl v s + y U2 v A J 

where v u is the vacuum expectation value of the higgs field h u . The other terms give the 
Majorana mass matrix: 




m M 



whose eigenvalues are: 



iu + 2bv s -bv s + cv A -bv s + cv A 
-bvs + cv A 2bvs + cv A au — bvs (14) 
-bvs + cv A au — bvs 2bvs + cv A 



Mi = au + 3bvs — cv A 

M 2 = au + 2cv A (15) 
M 3 = — au + 3 bvs + cv A . 

Using eq. ([9]) , we can derive the light neutrino mass matrix, diagonalized by tri-bimaximal 
mixing, whose eigenvalues are: 

m ( Vu Y ( 3 Vui Vs ~ Vu2 Va + yu * u "> 2 

\ A / au + 3bvs — cv A 
m 2 = _fM 2 (2y V2 v A + yu3 uY 
VA/ au + 2 cva 
fVu\ 2 (3jfo vs + Vu 2 va -Vu 3 u) 2 
\ A / au — 3bvs — cva 



m 3 = 



We see that the neutrino masses depend on six unrelated complex Yukawa parameters, 
which offer more freedom to tune mass differences and then recover the phenomenology 
associated to neutrino oscillation. Notice also that no sum rules can be found in this case 
among complex eigenvalues. 
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2.1 Charged leptons 

It could be easier to work in a basis where the charged lepton mass matrix is diagonal. 
In order to understand how this naturally arises in an S^-based model, we observe that a 
generic diagonal matrix mf (with diagonal entries different to each other) is left invariant 
under the action of an element A of S4 only if such an element is itself diagonal, with 
different phase factors at each diagonal entry, as it can be understood requiring that the 
relation 

A* mf f mf A = mf ] mf (17) 

is satisfied. The generator T, for example, is such an appropriate matrix. Since the 
charged lepton matrix is generated after spontaneous symmetry breaking, one could 
choose flavon fields with vevs invariant under the action of T. Instead, we prefer the 
choice 

(cp T ) = v T (0, 1, 0) (18) 

which, even breaking completely the group 5*4, not only guarantees the diagonal form of 
the mass matrix but also generates the hierarchy among lepton families without introduc- 
ing any addition U(1)fn symmetrjU. In the same way, we choose the vev for the doublet 
flavon 1] as: 

(V) = % (0, 1) ; (19) 
the corresponding lagrangian in the charged lepton sector is as follows: 

C = y -L T c {£vT)h d + 

^ c (£ (fr) h d + |f (V <Pt) h d + 

y^e c £[M^T)2}s 2 h d + |fe c h d + (20) 
|§e c £ [r, (WTkk h d + ^e c e[^ T ( W ) 2 ] 3a h d + 
|fe c £(A^ 5 ) 32 h d . 

It is interesting to observe that the last term in the lagrangian would be the dominant 
one and would drive the electron mass to a too large valued). However, the leading order 
structures of the vacua in eqs.f llOp - ffTT]) prevent this term to appear. We will show later 
that this is not the case when the next to leading order corrections to the flavon alignment 



2 This works in the same way as described in [5] noticing that, like in the case of A4, (0, 1, 0) 2 = (0, 0, 1) 
and (0,1, 0) 3 = (1,0,0). 

3 This is because this term explicitly breaks the residual symmetries needed to generate the correct 
hierarchies between the charged lepton masses. 



6 



are taken into account. After symmetry breaking, the mass matrix has the form: 



VVt I -&i [vT(Ve 1 + 2ye 2 ) -2v T v n y e3 + v*y e4 ] 

™>t = —f— i ( 2 Vi*i v t + y^ 2 v v ) o 

V Vr 

(21) 

where Vd = (hd). To estimate the order of magnitude of vt and v v , we can use the 
experimental informations on the ratios of lepton masses. Assuming that the combinations 
of the y coefficients are all of 0(1), one obtains: 

— ^ ] ~ 2e r + e„~0.06 



m T 



3e 2 T -2e T e v + el~ 0.0003 



where we introduced the small quantities 

e T = VT /A e v = v v /K. 
These relations are satisfied for: 

(kr|,|^|) ~ (0.017,0.029). (22) 

so that we can roughly assume that both et and e v are of the same order of magnitude, 
e ~ 0{\l). 

2.2 Superpotential and vacuum alignment 

The most general driving superpotential Wd invariant under S4 x Z 5 with R = 2 is given 
by 

Wd = gi((Po(p s <Ps) + #2(^o¥^K + #3¥?o(¥^A) + 

2 4 A (AA)+£ 5 Ao(<^< ? ) + <7 6 A A£+ (23) 

hi (pl((p T (p T ) + h 2 <Pq(t)<p t ) + n p (v ? TV ? r) + r 2 Po(w) 

= w L d o (^,A )+w L d °^l, P0 ) 

where all possible contractions among flavon fields are understood. The equations which 
fix the components of the vevs of the various flavon fields are obtained solving a system 
of equations obtained deriving Wd with respect to the component of the driving fields. 
The charge assignment reported in Tab{T] allows to separate the equations into two sets 
of independent relations among the fields appearing in the neutrino sector and in the 
charged one, respectively. For the latter we have: 
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dw d 
dw d 



-- 2h 1 {ip T \ - ¥t 2 V?t 3 ) + h 2 (Vi <Pt 2 + V2 <Pt 3 ) = 
= 2/ii(v9 T 2 - (fTi V?t 3 ) + h 2 (r]i +V2<^t 2 ) = (24) 

i 



2h 1 (<p T3 - (p Tl (p T2 ) + h 2 (vi ¥t 3 +r]2 ( PT 1 ) = 



'02 
9w 



dp 

whose solutions are: 



- ri (<£>t? + 2 v2t 2 V?t 3 ) + 2 r 2 r^i r? 2 



{//) = '•„ (0.1). = r r (0.1.0) . r„ = 2 ( ^ ) r, . (20) 



It is important to observe that the last equation is crucial to avoid another solution of 
the form (rj) = (1, —1), (<Pt) = (1, 1, 1)- We also recover the relation between the vevs vt 
and v v of eq. (l22l) . In the neutrino sector, the set of equations read as follows: 



dw d 

9<P S 01 
dw, 



dip s 



- 2g 1 ((ps 2 l - (fis 2 ¥s 3 ) +92^, <psi + 93 (Ai <Ps 2 + A 2 <psa) = 
2g 1 ((p s 2 2 - (p Sl (p S3 ) + g 2 ^cp S3 + g 3 (A 1 cp Sl + A 2 (p S2 ) =0 



y 02 

dw, 



dw, 



2 9i(<fisl- <Psi¥s 2 ) + 92^s 2 + 93(Ai¥s3 + &2¥si) = (26) 
- 5-4 + 5(5 (v? 53 + 2 <y?Sl V?5 2 ) + #6 A 2 £ = 



5A i 

= g± A2 + 05 (<^52 + 2 <ps! y?<? 3 ) + 5-6 Ai i = . 

C^02 

The system is solved by: 

(0=u, (A)=v A (l,l), W = (1,1,1), (27) 
with the additional relations 

v A = -— 28 
2 £3 

202 03 06 -02 04 \ 2 / on x 

7^ 2 M ( 29 ) 

1205 # 3 / 

and u undetermined. We have explicitly checked that the solutions of the vacuum align- 
ment equations are unique; in fact, requiring for the generic field $ LO to be shifted as 
$ LO + 5$, we found that the components 5<3? are all in the same directions of the corre- 
sponding $ LO ; thus, we do not need to introduce any soft term to drive the superportential 
into the wanted minimum. 
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3 Next to leading order 



The discussion of the corrections to the previous results starts with a study of the next to 
leading order structure of the vacuum alignments of the flavon fields. It will turn out that 
such corrections will be enough to guarantee deviation from TBM at a level compatible 
with the recent experimental results. 



3.1 Corrections to the vacuum alignment 

The next to leading order terms mix the charged lepton and neutrino sectors in a non- 
trivial way. We want to find perturbations of ( 12511271) of the form 

(ips) = (v s + Sv Sl ,v s + 5v S2 ,v s + 5v S3 ) 
((p T ) = (5v Tl , v T + 5v T2 , v T + 5v t 3 ) 

(A) = (v A + 5v Al ,v A + 6v A2 ) (30) 
(0 = u + 5u. 

On a general ground, we have eleven unknowns but only nine equations (six from the 
two triplets, two from the doublet and one from the singlet); then, we expect two of the 
previous shifts to remain unconstrained. To better understand the output of such an 
analysis, we study the (<^q , A ) and (pQ,p ) sectors separately. 



3.1.1 The (<#J,A ) sector 

This sector is responsible for the alignment of the fields (ps and A. At order 0(1/ A), 
the three-field terms entering the superpotential are combination of (p? and rj; collecting 
these terms, we have: 

5w d (<p%) = £? = i J <^0 (<PT<PT<fiT)i + S i=4^ ¥>0 + X ^° (fTfTV)i (31) 



5lU d (A ) = Ej=l-£ A (VVV)i + S i=3^ A (¥T¥TV)i + X A {¥T^T¥T)i (32) 

where, in both cases, the index % of the trilinear terms represents different S4 contractions. 
These NLO corrections have to be evaluated with the leading order vev's in eq. (l25|) . so 
that they do not contain any of the unknowns under investigation. In particular, the 
structure of the vacua (125]) produces 

5w d (A°) = 0; 

and the total LO+NLO part of the superpotential responsible for the alignment of ips 
and A is given by: 

w d {cp$, A ) = w%°{cp%, A ) + 5w d (cp$) . 
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Symmetry arguments allow to understand the structure of the solutions, whose detailed 
expressions can be found by explicitly solving the system of equations according to the 
procedure of Sect. fl2.2p . In fact, we see that, after symmetry breaking, the non- vanishing 
terms in 8wd(<fo) are proportional to the vector (1, 0, 0), which is left invariant by the S4 
element (TSTS 2 ); this has a 2 f> 3 symmetry which forces the new vev of (fs to have 
the form (vs + ^fs^fs + Svs 2 ,vs + 8vs 2 )- At the same time, the shifts 5v Ai should all 
be equal because of the vanishing NLO term 5wd{A )- An explicit computation confirms 
these speculations and we can write the new vacua in the following form: 



(p s = v s 



1 + B S (^e 
{l + B s (fYe J 



A 



l + A A (f) 2 £ 
l + A A (§r) 2 e 



(33) 



where the coefficients As, B$ and A A are linear combinations of leading and next to 
leading order coefficients and, for the sake of simplicity, we have introduced the small 
parameter e' = (<ps) /A ~ (A)/A ~ (£)/A- Notice that the three Svsi also depend on the 
undetermined parameter Su with identical coefficients, so that they can be readsorbed in 
the leading order result. 



3.1.2 The (y?o,p ) sector 

This sector is responsible for the alignment of the fields and rj. At order 0(1/ A), the 
new terms in the superpotential are combinations of (ps, A and £; we have: 



(34) 



t 



A 



A 



A 



and 



Sw d (po) 



^ Po ivsVsVs) + ^Po (A<Ps<Ps) + ^Po (Vs <Ps) £ + 
^Po(AAA) + ^p (AAK + ^Po£ 3 . 



(35) 



In this case, the corrections to <pg do not conserve any of the S4 basis elements and we 
expect that the components of (<Pt) wm point to different directions. Also the (77) 's vev is 
not preserved at the next to leading order and, consequently, the shifts 5v r/i are different 
from zero and different to each other. We choose to treat Sv m as the second undetermined 
parameter, so that the new vacua can be cast in the following form: 



/ 



+ Br Ci) 



e 



\ 




(36) 
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Notice that a correction proportional to 5v V2 also appears in the second component of ip? 
in the form ip T2 = vt — (h 2 /2hi) 5v V2 + 0(e') which, using the last relation in eq. fl25|) . 
gives ifT 2 = — (h 2 /2h 1 ) (vt + Sv m ) + 0(e'), so that the shift Sv m can be readsorbed into 
a redefinition of vt- 

3.2 Charged lepton mass matrix 

The next to leading order corrections to the mass matrix of eq. (!2ip come from the cor- 
rections to the vacuum alignments in eqs. (l33l) and (136|) and from higher dimensional 
operators, suppressed by a relative 0(1/ A) with respect to each of the terms quoted in 
eq. fl20l) . In the following we will study them separately. 

3.2.1 Corrections from vacuum alignment 

The main features of these corrections are related to the fact that all the vanishing entries 
in the matrix of eq. (121 j) are suppressed by one additional Oil/ A) factor compared to the 
diagonal entries, except for the first line, where the 0(1/ A 2 ) operator in eq. fl20|) is non 
vanishing. As a result, we get the following charged lepton mass matrix at the NLO: 

(3 '2 '2 \ 

aiE A a 2 ee —a 2 ee \ 

biej b 2 £ 2 b 3 e' 3 (37) 

ci ^ c 2 ^ c 3 e ) 

where the coefficients a iy bi and c, can be easily reconstructed from eqs. (l2"Uj) and ( 13T)j) . The 
matrix m\ mi can be diagonalized by the unitary transformation 



Us 



( 1 (h, (£2, dlv \ 

1 U 2 e 2 > \ c 3 e 2 I 

_6i e 3 y (£2. £ 3 )* 

\ C3 E 2 C3 e 2 / 



(38) 



so that, at the lowest order in e', the coefficients of the electron row in eq. (|371) do not 
contribute. From the matrix in eq.f l38|) . we can compute the £7pmns mixing matrix from 
the lepton sector only: 

f^PMNS = U\ UtbM — 

( V / i + 75 ^ + ^ (£)* 73-73 ^ + ^ (£)* 75 ^ - & *) (£ 

-yfl + y/l* \/i + 73 a (#)-73 c *(#r -y/l 



where a = bi/b 2 , b = Ci/c 3 and c = c 2 /c 3 . Then, any entries of £/tbm get corrected by 
terms of 0(e' 3 /e 2 ) which can be as large as 0(\c) to fit the experimental data. 
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3.2.2 Corrections from higher dimensional operators 

For the electron case, there are many of such operators, coming from the following con- 
tractions: 



(AVt) (A 3 r?) (<p 3 g(p T ) [(p 3 s ri) (A 2 W5 ) 



(AVO (AVst?) (A 2 ^) (^ t A) (<^t£) (40) 

(f 2 s Ar]) (£, 2 (p T A) (£, 2 ipsv) {£, 2 <Pt<Ps) 

It is easy to understand that they contribute to the electron row by terms of order (9(1/A 4 ), 
which are unimportant because they do not sizebly modify the diagonalizing matrix Ug. 
Also, the corrections to the r row are of 0(1/ A 3 ), that is the coefficients b and c appearing 
in the Upmns are modified at the next to next to leading order (NNLO) and can be safely 
neglected. For the muon case, the following higher dimensional operators modify the 
coefficient a appearing in eq. fl39|) at the NLO : 



_ V.. NLO t V., NLO 

C NL0 = Zl 1 ^» c £(tp sW s)h d + Zt = 3-^ f , c e(<psAA)h d + (41) 
Vj ^^Vse+ Vj ^^^sA)i. 

In conclusion, both types of corrections (from vacuum alignment and from higher dimen- 
sional operators) to the charged lepton mass matrix contribute to generate deviation from 
TBM at O(Xq). 



3.3 Neutrinos 

Also for neutrinos, we have to take into account corrections from both vacuum alignment 
and higher dimensional operators. In particular, the NLO corrections of the vevs of the 
flavon fields affect both Majorana and Dirac masses with the same pattern, as we can see 
from eq. ()12p . Different higher order corrections are generated by higher order operators. 
These corrections turn out to be negligible compared with the NLO and of the same order 
of magnitude as those induced by the Weinberg operator (see later). 



3.3.1 Corrections to the Majorana and Dirac mass matrices 

The relevant contributions come from the NLO vacuum alignments. In particular, eq.f lT4|) 
is modified only by (<ps) at the NLO because the other two vevs, namely those of A and 
£ entering in the lagrangian in eq. (fT2|) . are aligned along the LO direction. However, 
the shifts in (ip s ) are not enough to introduce independent corrections to all the (six) 
elements of the Majorana and Dirac mass matrices; in fact, the explicit structure of the 
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NLO corrections are as follows: 



£' 



Sm M = b - ev s S (42) 
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5m D = v u y—j e y Ul S (43) 



where the common matrix S is: 

— Bs — Bg 

5=1 -Bg 2Bg -Ag J . ( I I) 

-A s 2B S 

With these results, we can build the light neutrino mass matrix; since the NLO matrix 
elements have now complicated expressions in terms of the leading order Yakawas of 
eq. (fl~2l and of the correction coefficients As and Bs, we prefer to summarize its structure 
as follows: 




X 


y 


y \ 




f x' 


y' 


y' \ 


y 


X + V 


y-v 


HI 


y' 


z 




y 


y-v 


X + V J 




U 


z' 





might = £ \yx + vy-v\+e\y'zz'\, (45) 



where the relation z+z' ^ x'+y' implies that m Ught is not diagonalized by Utbm- However, 
it is easy to show that the almost symmetric structure of the e 3 contribution prevents 
to generate corrections to the TBM values of 9 is and #23 because the diagonalization of 
mught is achieved by a matrix U v = Utbm + SU, where SU has a vanishing last column. If 
we want to obtain a non-vanishing # 13 , we need to distinguish the elements (12) from (13) 
and/or (22) from (33) in eq. f)44p . Without taking into account the corrections from the 
charged lepton sector, this can be accomplished only at the NNLO. In fact, a Z 5 singlet 
built with the v c v c bilinear requires a charge u for the accompaining fields, and this can 
be achieved with at least three fields, which induce terms of 0(1/ A 2 ). There are more 
than twenty of these terms, generated by the following contractions: 

(AV) (A 2 r?) (vfapr) (cp 2 sV ) (£V) (fr,) 

((p T A(p s ) (<p T AO (<Pt<Ps€) (^Vfs) (A77O {tpsvO ■ 

Similarly, the Dirac mass matrix receives corrections at the same NNLO relative to the 
leading terms due to the fact that the term iv c €) h u has a total Z5 = u 4 . 

Beside the previous operators, we also have to consider effective terms of the form 
llh u h u ; given its charge assignment (Z 5 = a; 4 ), the Weinberg operator arises at 0(1/ A 2 ) 
with the insertion of one flavon field, as for the Dirac and Majorana terms: 



Wf = £ (ai<f S + a 2 A + a 3 £) (t h u t h u ) . (46) 
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After spontaneous symmetry breaking, W£ s generates terms which are of order: 



(47) 



Compared to the NLO corrections of the Dirac mass term, eq. (j43l . we see that: 

/r U «M 1 

e 

then this type of operators is more suppressed with respect to the Dirac mass corrections 
(and even more if compared with the Majorana mass terms) and can be safely neglected. 
In conclusion, from the neutrino sector only, the following mixing matrix arises: 



2 (x'+y'-z-z')* 
3_ 9y* 

2 (x'+y'—z—z')* 

3 9y~* 

2 (x'+y'-z-z 1 )* 

3 9y~* 



j_ , 2 (x'+y'-z-z') I ; ' 

(x'+y'-z-z') 

9V3y 

(x'+y'-z-z') 
9V3y 



^) o 



\ 



2 



(48) 



3.3.2 Mixing angles at the NLO 

From the previous discussions, it clearly appears that the NLO corrections to #13 and #23 
come only from Ug whereas 812 is modified by both Ue and U v . They are as follows: 



S13 



S12 



S23 



\U, 



<>3 1 



1 /e 3 

T2 {A - BY (f 



It/, 



r;2 



4= -^(^ + 5)' 



-4 + 1 



v/2 



where we used the short-hand notation 



2 (at + y'-z- z') 

3 % 



(49) 



In particular, the first relation can be used to put a bound on |e'|; in fact, given the 
maximum allowed value for #13 and assuming that \A — B\ ~ 0{1) we get: 



< 



-,1/3 



(50) 



which is close to (2 A^) for 6f^ ~ 10°. 
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4 A bit of phenomenology 



The model we have presented has a huge parameter space, made by the six complex 
Yukawa couplings y Vi and a, b and c. It is clear, and we have checked this numerically, 
that the various phenomenological constraints, coming for example from the smallness 
of the parameter r = Am^ oZ /| Am^ im |, are easily satisfied, for both type on neutrino 
mass hierarchies. It seems then more interesting to ask whether the model can still 
give an acceptable phenomenology in some particular cases, like for example allowing all 
the Yukawa couplings to be of 0(1) and the flavon vevs (<ps), (A) and (£) to be of the 
same order of magnitude^. In this case, we do not expect any huge hierarchies among 
the heavy neutrinos, see eq. (fT5]) . Our numerical study aims to predict some interesting 
physical quantities only imposing the following 3-<r experimental constraints 



| Am 2 
r 



atm I 
exp 



> 



2.41 ±0.34 x 10 
0.032 ±0.006 



-3 



eV 2 



(51) 



also taking 



rrii 



< 0.5 eV. 



their sum is shown in Fig.( 



The resulting spectrum of the light neutrino masses and 
]). In the left panel, we chose to plot the ratio |m3|/|m2| 



s I s 
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normal hierarchy **"»-». 
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io- 



Figure 1: Neutrino mass spectrum and sum of neutrino masses, as predicted from the 
model with all Yukawa couplings of 0(1) and the flavon vevs of the same order of mag- 
nitude. Left panel: behaviour of the ratio \m3\/\m2\ as a function of the lightest neutrino 
mass \mu g htest\, for both neutrino hierarchies. Right panel: sum of the light neutrino 
masses as a function of \mughtest\- Also shown are the bounds from fTbJ (upper solid 
lines) and from [Tfif + J7?[ / (lower solid lines). The vertical line in both panels is the 
future sensitivity of 0.2 eV on mughtest from the KATRIN experiment [Tffl . 

as a function of the lightest neutrino mass \Tnu g htest\, which is |mi|(|ra3|) for the normal 
(inverted) hierarchy. We clearly see that the largest hierarchies, obtained at the smallest 

In the numerical simulation, this means that we allow the absolute values of the Yukawas to be in 
the interval [1/2, 3/2], while no restriction whatsoever has been imposed on the flavon vevs a part from 
being all equal within a factor of 10; also the large scale A is left free. 
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allowed \mu g htest\i are at the level of O(10) for the normal ordering and O(10~ 2 ) for the 
inverted one. The ratio tends to a degenerate spectrum in both cases for \m Ughtest \ ~ 1CT 1 
eV which, however, is disfavoured in our model. In the right panel we show the sum of the 
light neutrino masses as a function of the lightest neutrino mass mu g htest- The vertical line 
denotes the future sensitivity of 0.2 eV on of \mu g htest \ from the KATRIN experiment [Hj, 
and the horizontal lines are the cosmological bounds [15] at 0.60 eV, obtained combining 
the data from ref. [16], and at 0.19 eV, corresponding to all the previous data combined to 
the small scale primordial spectrum from Lyman-alpha (Lya) forest clouds [T7j . We see 
that our model predicts £ m ; too similar for both hierarchies to be distinguished using 
the current cosmological information on the sum of the neutrino masses; however, such a 
discrimination could be possible if some improvements on these bounds would be achieved 
in the near future. 

Finally, we present in Fig.(T5]) the predictions for the values of the effective mass \m ee \ 
as a function of the lightest neutrino mass, for both normal and inverted hierarchy. We 
also show the future sensitivity of the KATRIN (vertical solid line) and of CUORE [IB] 
(horizontal solid line at 15 meV) experiments. The main feature of the analysis is that 



1 



0.1 



— 0.01 



0.001 



10- 



inverted hierarchy 



normal hierarchy 



0.001 



0.01 0.1 

\mu g htest\ (eV) 



CUORE 



Figure 2: \m ee \ as a function of the lightest neutrino mass, for both normal and inverted 
hierarchy. The filled regions correspond to the possible values of \m ee \ in the limit of exact 
tri-bimaximal mixing, with mass differences computed at the central values in eg. 137]) . The 
horizontal and vertical lines are the future expected bounds on \m ee \ and \mu g htest\ from 
the CUORE and KATRIN experiments, respectively. 

a large set of points falls into the region of \mu g htest\ around 10~ 2 eV but, given the still 
large number of parameters of the model, many values of \m ee \ can be obtained for both 
hierarchies among the experimental allowed ranges. On the other hand, the region above 
this value (the regime of degenerate spectra) is strongly disfavoured. 
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4.1 Leptogenesis 



The formal description of the asymmetry parameters can be done in the general con- 
text of arbitrary Dirac mass matrix because the final expressions are quite compact and 
transparent. The asymmetry parameters are defined as follows: 



1 



i7r(yyt) i 



Im 



ij 



f 



I M« 



(52) 



where the hat matrices are Yukawa matrices evaluated in the basis in which the Majorana 
mass matrix is diagonal and Mj are the Majorana masses. For supersymmetric theories, 
the /-function is given by: 



/(*) 



x 



2 , fl + x 
+ log 



X 



(53) 



Defining Q as the unitary matrix which diagonalizes the Majorana mass matrix, the LO 
Yukawa matrix in this basis is given by 

v u Y = v u f2 T Y v = il T m d 

and the product YY^ reads: 

YY j = tt T Y v Yl tt* . (54) 

At LO, Q = Utbm and the product YY^ is a diagonal matrix: the parameters are all 
vanishing [12]. At the next to leading order, one has to take into account the corrections 
to the Yukawa matrix as well as to the Majorana mass matrix, which reflects in a different 
structure of both the Q and Y v matrices in such a way that: 



v u Y v = m D 



e 2 y Vl S) =v u (Y LO + 6Y) 



(55) 



where 5Q has a structure similar to the corrections 5U computed in Sect. (13.3. ip and 5Y 
is of 0(\q) compared to its leading order result. This means that the correction to the 
matrix product YY^ is given by: 

<f(YYt) = (5n) T Y LO Y ] LO n* + n T YjjoYIq (5ft)* + ft T S(YX) ft* • 

The first two terms do not contribute to the q parameters because they are complex 
conjugate of each other and do not contain any imaginary part; then the only contribution 
arises from the last term. In the basis in which the charged leptons are diagonal and 
considering e as a real variable for simplicity, one easily obtains: 



(56) 




V2e- i( -^- 
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where the a, coefficients are complicated functions of y Ui , As and B$- Then, at leading 
order, the e parameters are given by: 



n = 1^1 - 



2 4 (a s - B s f \.yvi(y»i + + 3 l^i I 2 ] 2 / (pfp) sin [2(0i - 02)] 



e 



2 



(A s - B s f \y^^ + 2 ^ 3 ) + 3|^J 2 ] 2 / (pg$J sin [2(0i - 2 )] 



U^'^tt^ 4 |^ 2 |2 + |^ 3 |2 + V 4(Mi/3) ( 57 ) 

e 3 = 0. 

A relevant feature of the model is that e 3 always vanishes; however, we can see that, 
barring possible fine-tunings in the parameters and/or suppressions or enhancements due 
to /(|a^| 2 ) sin [2A0], ei^ can be of the right order of magnitude to fulfill the experimental 
requirements for a successful leptogenesis because 

ei,2 ~ Ac ~ 6 x 1CT 6 . (58) 



5 The quark sector 

It is well known that the extension of a flavour symmetry from the neutrino sector to 
the quark one is highly non trivial due to the "non-trigonometric" structure of the quark 
mixing matrix Vckm- The common way to introduce quarks in these kind of models 
is to put the two SU{2) doublets of left-handed quarks into a doublet representation of 
S4 and assign the heaviest ones (top and bottom quarks) into singlets, in such a way 
to easily maintain the hierarchy among mass eigenstates. Better results are obtained if 
other extra symmetries, i.e. Froggatt-Nielsen or extra Zn, are introduced in order to 
further suppress the unwanted couplings; this is particularly true for the m up /mt p mass 
ratio, which otherwise tends to be larger than the experimental counterpart. It is clear 
that such scenarios are more flexible than the simple S4 x Z 5 illustrated in this paper; 
however, it is important to stress that it is still possible to get a satisfactory description 
of the quark sector without invoking any other extra symmetries and using the triplet 
representation of S4 instead of the doublet ones. In the down-quark sector the simplest 
choice is to copy the coupling allowed in the charged lepton sector because the hierarchy 
between the d, s and b quark masses is quite similar to that existing for e, /x and r 
leptons. On the other hand, the mass hierarchy in the up-quark sector does not follow 
such a simple prescription and we can arrange the S4 x charge assignment in several 
ways. One possibility is summarized in Tab.(j2J), where Q is a triplet of SU(2) left-handed 
doublets. The corresponding leading order lagrangian is the following: 
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Field 


Q 


d c 


s c 


b c 


u c 




t c 




3i 


1 2 


ll 


ll 


h 


ll 


ll 




1 


u 3 




u; 




1 


a; 4 


U{1)r 


1 


1 


1 


1 


1 


1 


1 



Table 2: Transformation properties of quarks under S4 x Z 5 and £7(1 )_r 



Vb 
A 



. 2 s c (g <^ T <p T ) h d + || s c g (77 ip T ) h d + 

^|cf Q [<Pt(¥t¥t)2}32 h d + ^jd c Q [^(^Wrkk h d + 

^|cf Q [t) {<p T <p T ) 3l } 32 h d + ^|d c Q [<p T (7777)2)32 + 

||d c Q (A^ 5 ) 32 ^ + |t c (g^ T ) ^ + (59) 

||c c (g Vt) K + V -^c c Q (77 ^t) K + 

-^u c Q[A(^ T ^ T ) 3l } 3 , 2 h u + j§u c Q[^ s (^tVt)2}3 2 K + 

j^u c Q [ip s (p T ^ T ) 3l } 32 h u + ^-m c g [ip s (w)2h 2 K + 

j^u c Q [77 (y? T y? 5 ) 3l ] 32 /i u + j^u c Q [77 (vWskk /i u . 

This lagrangian gives rise to a diagonal mass matrix for the down-type quarks and an 
almost diagonal mass matrix for the up-type quarks (all the up entries are filled with terms 
of the same order e 2 e'). This picture slightly changes when including the next to leading 
order effects. As for the case of charged leptons, in the down sector the relevant NLO 
corrections come from the corrections to the vacuum alignment of the flavon fields, which 
are all of relative 0(1/ A) with respect to their leading order counterparts, whereas the 
higher order operators give corrections of the same size only for the s-quark, via couplings 
like 



(tp 3 s ), (y? s A 2 ), (^ 2 ), (PsAO 



(60) 



and are much more suppressed for the down and bottom quarks. Similarly, for the up- 
type quarks the mass matrix is modified by the NLO structure of vacuum alignment of 
the flavon fields and, for the c-quark, also by the same higher order operators modifying 
the s-quark entries in eq.( |60l . All in all, the mass matrices for both type of quarks are 
given by: 



m do 



afe 3 
bfe' 3 

„d e' 3 



d '1 

a%ee A 



,d n 2 
'3 



b^e 



„d £ 



d '2 

— a% ee z 
b d 3 e' 3 



c, e 





( a\e 2 e' 


al e 2 e' 


a 3 l e 2 e' 




b^ e' 3 


b u 2 e 2 


b\ e' 3 






'3 

C 2 e 





(61) 
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with the following mass eigenvalues: 



m d 
m s 



v d 



afe 3 + 



a 2 e e' 3 



v d b d 2 e z + 0(s 4 
v d 4e + 0{e 4 ) 



(62) 



and 



m c 
m t 



a u, £l+ (bjal_a 1 c l \ 
V b\ c\ ) 



v u b u 2 e 2 + 0(e' i ) 
v u c u 3 e + 0(e' 4 ) 



(63) 



The previous mass matrices can be diagonalized by unitary matrices is such a way that: 



m\ m d ) diag 



where: 



b± £f 
b| e 2 
4 e_ 



r d '3 



4 e± 
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4e± 

4 e2 



U j UL (ml m u ) U UL 

U d r ( m d m d)U dL , 



t 



u, 
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1 

bj_ {_ 
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c l £ 





(of 


'3 


^ J 


\4 


e 2 


1 
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\4 


'3 

E d 
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_£i il 

n u _2 
c 3 £ 




1 



(64) 



The resulting Vckm is given by: 



VcKM = Ul U dL 



4 




e d 


4 




£ 2 


L l 




£ d 


c 3 




£ 2 



\ 



As for other flavour models, the matching of the Vckm and the quark mass ratios to their 
experimental values requires some fine-tunings between the Yukawas. As anticipated, all 
the experimental mass ratios in the down sector are easily reproduced for the natural 
values af, b 2 , c d ~ 0(1) because m d /m s ~ m s /mf, ~ e ~ \ 2 C . A moderate hierarchy is also 
present in the up sector and it mainly depends on the parameter e', which we estimated 
in Sect. (13.3.21) to be not larger than some units of A^. This same parameter appears in 
the off-diagonal entries of Vckm with the same power in each entries. This is essentially 
the reason why it is difficult to explain at the same time the mass hierarchy in the up-type 
quark sector and the off-diagonal values of the quark mixing matrix. For example, one 
can treat e' and tan j3 as free parameters of the model. In that case, one can use tha ratio 
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m u /m c to fix the value of e' and m^/m^ for that of tan/3; assuming for all Yukawas the 
natural 0(1) value, we get: 



m c 



m d 



tan j3 



and then 



— - ) tan (3 = X 

"Wexp 



12. 



(66) 



(67) 



cxp 



It is easy to verify that the previous results are enough to accomodate all the independent 
mass ratios that can be built from six different quarks. However, e' turns out to be very 
small and, having assumed all Yukawas of 0(1), it is really difficult to enhance the off- 
diagonal elements of Vckm to their experimental values. The other possibility is to 
preserve e' ~ 0(X C ), requiring that 



6f 

hi 



b\ 



0(1/X 
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\(4 
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0(X 
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sy 


\4 





0(1). (68) 



This fine-tuning is a condition on the individual Yukawas, which can be chosen to be: 



c 2 



bf 



b\ ~ 0{1) 
0{Xc) 
0{l/X c ) , 



(69) 



providing that b% , ^> 1 to fit the charm and top quark masses. Although this situation 
is not completely satisfactory, it illustrates a way to account for many experimental infor- 
mations in the quark sector using a simple and unconstrained S4 x Z 5 model, the prize to 
be paid being a fine-tuning in four of the Yukawa couplings appearing in the lagrangian. 

Notice that it is not more difficult to modify the charge assignment proposed in Tab. ([2]) 
in such a way to obtain at LO a Vckm matrix with entries of 0(1) in the Cabibbo 2x2 
submatrix. This can be accomplished slightly changing only the Z5 charge of the quark c. 
In fact, giving the assignment c ~ (li,u; 3 ) under (S4, Z 5 ), the c-quark part of lagrangian 
is now: 



||c c (Q cp T A) h u + ||c c Q (<p T <p s ) h u + ||c c Q <ptZK + ^c c Q {<p s rj) h u (70) 



and the c-quark entries in the mass matrix are modified according to 
(c c ,w) :&?£■£■' (c c ,c):b%ee' (c c ,t) : b% e e' . 
Correspondingly, the u and c-quark masses are: 

b\ a u 2 



Oi — 



1M 
°2 



e 2 e' 



v„ bo e e' 



(71) 



(72) 
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and the Vqkm has the following structure: 
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J 



(73) 



where, for simplicity, we introduced the short-hand notation K = a/ (&" ) 2 + (&f) 2 an( i 
considered real vevs and Yukawa couplings. In this case, one can reproduce at the same 
time the correct values for the (11) and (22) entries of V C km as well as the mass ratio 
m u /m c imposing b% ~ O(Xq); playing with e' and the other Yukawas one can also repro- 
duce the off-diagonal entries, but one or more fine-tunings in the coefficients in front of 
e 3 /e 2 ratio are obviously needed. 



6 Conclusions 

We have presented and discussed an £4 model for TB mixing (of the see-saw type) and 
quark mixing which, in spite of being based on a most economical flavour symmetry and 
field content, it is still phenomenologically viable. In the neutrino sector, we realized the 
most general neutrino mass matrix diagonalized at LO by TB mixing. At the NLO, all 
mixing angles receive corrections at the level of O(Xq) and, in particular, the elusive #13 is 
predicted to be within the sensitivity of the experiments which are now in preparation and 
will take data in the near future. Also, we have shown that the leptogenesis parameters 
61,2 can be of the right order of magnitude to satisfy the requirements needed to reproduce 
the observed asymmetry and that an acceptable phenomenology is still obtainable even 
with all Yukawa couplings at 0(1). 

In the charged lepton sector, the mass hierarchy is determined by the S4 x Z 5 flavour 
symmetry itself without invoking a Froggatt-Nielsen U(l) symmetry. The NLO correc- 
tions to the mass matrix turn out to be relevant for a #13 7^ 0. 

In the quark sector, we have discussed in detail to which extent the model can repro- 
duce the data on masses and mixing, emphasizing where it is successful and the reasons 
for its failures. Two different LO examples for Vckm have been given, one proportional to 
the identity matrix and the other with 0(1) elements in the 2x2 sector; both realizations 
need almost the same amount of fine-tuning to reproduce the off-diagonal entries at the 
NLO. 
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A The Group S A 

We adopt the following convention for the generators S and T, according to [9] 



2\2 



S 4 = T 3 = (ST 



In the different representations, they can be written as reported in Tab. ([3]): 



(74) 
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2lu 2cu 2 -1 
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Table 3: Generators S and T in different representations . 
The 24 elements of the group belong to five conjugacy classes 
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The explicit expression of the elements in the 2-dimensional representation is: 
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while for the 3-dimensional representation 3i the elements are 
• Ci : I = 
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2co 2co 2 
• 2 2 



-CO 

2 



-U! 



T 2 S= - 



-1 2co 2co 2 
2co 2 2 -co 
2lo -co 2 2 



1 
7',S7',S' 2 = | 1 
1 



STS 2 




S 2 TS = 



• C, : 



S = - I 2co 2uo 2 
3 1 o . ,2 



ST= - 
3 




2^ 





/ -1 2w 2 2w 
2co 2 2co -1 
\ 2w -1 2w 2 




S 3 T 3 



-1 2co 2co 2 

2 2co -co 2 

2 -co 2co 2 

-12 2 

2co 2 2co 2 -co 2 

2uo -co 2uo 
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For the 3-dimensional representation 32, the matrices representing the elements of the 
group can be obtained from the list for the representation 3i in the following way: for 
Ci,2,3 are the same, while for C 45 are the opposite. 

In the previous basis, the Clebsch-Gordan coefficients are as follows (a* indicates the 
elements of the first representation of the product and j3i the second one): 



li <S> rj 



7] <S> li = rj with rj any representation 



1 2 
1 2 



li ~ a : 
2 - 



a Pi 
-a(3 2 



1 2 
2 

3i 



I2 <S> 3 2 = 3i ^ 

The multiplication rules with the 2-dimensional representation are the following: 

li ~ a 1 f3 2 + 0! 2 /3i 




2©2 = liffil 2 ©2 with < 



1 2 ~ aifc - a 2 fti 



2 © 3i = 3i © 3 2 



3i 



with < 



2 © 3 2 = 3i © 3 2 



3i 



with < 



aciP 2 + Oi 2 f3 3 
ai/?3 + (X2P1 
eeiPi + a 2 (3 2 

a>i/3 2 - a 2 /3 3 
ai/3 3 - a 2 f3i 
<*iPi - o: 2 /3 2 

a\(3 2 — a 2 (3 3 
ai/3 3 - a 2 /3i 
aiPi - a 2 (3 2 

ui/3 2 + a 2 (3 3 
(X1P3 + <*ifi\ 
aiPi + ® 2 (3 2 
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The multiplication rules with the 3- dimensional representations are the following: 
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